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The intriguing appeal of circuits lies in their modularity and ease of fabrication. Based on a toolbox of simple 
building blocks, circuits present a powerful framework for achieving new functionality by combining circuit 
elements into larger networks. It is an open question to what degree modularity also holds for quantum circuits 
- circuits made of superconducting material, in which electric voltages and currents are governed by the laws 
of quantum physics. If realizable, quantum coherence in larger circuit networks has great potential for advances 
in quantum information processing including topological protection from decoherence. Here, we present theory 
suitable for quantitative modeling of such large circuits and discuss its application to the fiuxonium device. 
Our approach makes use of approximate symmetries exhibited by the circuit, and enables us to obtain new 
predictions for the energy spectrum of the fiuxonium device which can be tested with current experimental 
technology. 



In the search for a viable architecture for solid-state quan- 
tum information processing, superconducting circuits have 
been the focus of immense interest™. While research efforts 
have led to the remarkable improvement of coherence times 
by nearly 5 orders of magnitudes 5 relative to those in the pio- 
neering experiments a decade agcPEl superconducting circuits 
have remained extremely simple - especially when compared 
to circuits found in commonplace electronic devices. Whether 
phase, flux, or charge qubits, most superconducting circuits 
consist of less than a handful of circuit elements. 

Experiments with the fiuxonium device - a superconduct- 
ing circuit with more than 40 Josephson junctions - have 
shown that a larger number of Josephson junctions, and hence 
degrees of freedom, is not necessarily penalized by reduced 
coherence times 89 . Experimental studies of linear Joseph- 
son junction arrays have advanced at a rapid pace^^l How- 
ever, despite considerable theoretical wortfGKD methods for 
detailed modeling of larger circuit networks are needed to suc- 
cessfully chart the future territory of quantum coherence in 
networks of increasing size to, e.g., further explore the pos- 
sibility of topological protection from decoherenc o 28 '' 3 1 1 The 
description presents a considerable challenge to theory due 
to the combination of several factors: the non-linearity in- 
duced by Josephson junctions, the increased number of low 
energy degrees of freedom, and the peculiar interactions be- 
tween them induced by flux quantization. As a key step for 
mastering these difficulties, we here present theory for the 
fiuxonium device starting from its complete degrees of free- 
dom. We demonstrate that circuit symmetries play a crucial 
role in the organization of the excitation spectrum and, em- 
ploying the relevant collective modes and their approximate 
decoupling 26 , we obtain a systematic approximation scheme. 

Non-linearity, interactions and a large number of degrees of 
freedom are challenges commonly encountered in the study 
of many-electron atoms. Our symmetry based approach re- 
sembles methods familiar from atomic and molecular physics 
where the weak breaking of symmetries leads to the well- 
known lifting of degeneracies in the fine and hyperfine struc- 
ture of spectra 3 ^! For the fiuxonium circuit, we demonstrate 
that approximate symmetry under the unitary group and un- 
der permutations of junction variables divide the excitation 
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FIG. 1: Circuit of the Fiuxonium device: Josephson junction array of 
nominally identical junctions with Josephson energy Ef and charg- 
ing energy Eq, shunting one weaker (Josephson and charging energy 
i5j b and Eq, respectively. 



spectrum into nearly degenerate subspaces. For realistic pa- 
rameters, the careful study of perturbations allows us to refine 
our description and provide new predictions for the collective 
excitations of the circuit. 



I. FLUXONIUM CIRCUIT 

The fiuxonium device 8 (Fig. |T|> consists of a Josephson junc- 
tion array with a large number N 3> 1 of nominally identi- 
cal tunneling junctions. One additional smaller junction (the 
"black sheep") shunts the array. The superconducting loop 
formed this way can be biased with an external magnetic flux 
$ ext , making the energy spectrum tunable. As typical of su- 
perconducting circuits, the nature of eigenstates and their de- 
tailed energy spectrum are governed by the competition be- 
tween charge transfer across junctions due to Cooper pair tun- 
neling and charging effects due to excess electric charge on 
individual islands. While the former favors charge derealiza- 
tion and definite phase differences across each junction, the 
latter promotes charge localization with definite Cooper pair 
numbers on each superconducting island. The ratio of Joseph- 
son energy (Ej) to charging energy (Ec) of the involved junc- 
tions quantifies this interplay. We use superscripts "a" and "b" 
in the following to distinguish between array (Ef /E c 1) 
and black-sheep parameters (E} ~ E^). 

The energy spectrum and corresponding eigenstates of the 
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superconducting circuit C are governed by the stationary 
Schrodinger equation He | ip ) = E \ ifi ), in which the circuit 
Hamiltonian is obtained from the Lagrangian C = T — U 
by circuit quantizatio n 33 ! 34 ! For each junction, Josephson tun- 
neling produces a potential energy term U m — —Ej m cos 9 m 
where 8 m denotes the phase difference across junction m. The 
dominant kinetic energy contributions arise from the charging 
of junction capacitances, T m = \G} m V^ m , where the voltage 
drop across junction m is linked to 9 m via Josephson's phase 
evolution equation 9 m = 27rVj„ l /$o and charging energies 
are related to capacitances via Eq = e 2 /2C. 

Fluxoid quantization^ 3 -^ forces the total magnetic flux 
through the loop to be quantized in units of the superconduct- 
ing flux quantum $o. This leads to the constraint X)m=o @ m + 
(p exl — 2irz where ip ext = 27r<J > e xt/ < I ) o is the phase offset due 
to external magnetic flux and z is an integer. The constraint 
reduces the number of independent coordinates by one and in- 
duces coupling among the remaining degrees of freedom. To 
incorporate the constraint while maintaining symmetry among 
array junctions, we eliminate the black-sheep variable 6> and 
obtain the Lagrangian 
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where, as a convention, sums over Latin indices always run 
over the range 1,...,N. The capacitive term involving the 
matrix Q describes the effects from capacitances between su- 
perconducting islands and ground (see Supplementary Mate- 
rial). 

To illustrate the content of equation ([TJ, it is instructive to 
note that C describes a single fictitious particle inside a peri- 
odic potential, albeit in A-dimensional space with N 1. 
Alternatively, it can be interpreted as a description of N dis- 
tinguishable particles, each moving in a Id periodic potential 
but subject to a peculiar interaction of collective type induced 
by flux quantization. 

The central idea of our approach in the following is to 
harness the large amount of symmetry present in the domi- 
nant terms of equation ([T} 43 . In particular, if ground capac- 
itances are negligible and if all array junctions possess the 
same charging energy E^ and Josephson energy E^ = Ej 
then L is manifestly Sat symmetric: any permutation a E Sn 
of the array variables leaves the Lagrangian invariant for any 
value of the external flux. We will refer to this idealization as 
the Symmetric Fluxonium Model. 

In non-relativistic quantum mechanics, such discrete sym- 
metries generally lead to degeneracies which are governed by 
the irreducible representations of the symmetry group. The 
simplest irreducible representations of the symmetric group 
Sat are the trivial and alternating representations familiar 
from particle and many-body physics where they dictate the 
symmetry of wavefunctions for indistinguishable bosons and 
fermions. In the case of superconducting circuits, degrees of 
freedom referring to different junctions generally remain dis- 
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FIG. 2: Normal modes for 
the Sjv symmetric fluxonium 
model. Plots show the array 
variable amplitudes 8 m for 
each normal mode, (a) Su- 
perinductance mode [(f)], for 
which all array junction am- 
plitudes 9 m are identical, 
(b) Difference modes 
for all of which the amplitude 
sum exactly vanishes. 



tinguishable, and the full plethora of irreducible representa- 
tions of Sat is realized. In this sense, the SFM constitutes an 
intriguing example of a many-body system with degenerate 
eigenstates that obey novel permutation symmetries beyond 
those of bosons and fermions. 



H. Sat SYMMETRIC FLUXONIUM MODEL 

From the circuit Lagrangian ([T} we now extract the rele- 
vant collective modes^ governing the low-energy physics, 
and discuss their connection with effective models employed 
in previous worlP 3 -^. A key ingredient in the construction 
of the low-energy modes is the observation that array junc- 
tions in fluxonium are dominated by Josephson tunneling, 
E^/Ef -C 1, while the black-sheep parameters and Ej 
are both roughly of the same order as the array charging en- 
ergy. 

For large arrays with junction number A ^> Ej/Ej, the 
potential energy U exhibits deep minima at positions where 
all array coordinates have the identical value 
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with integer z satisfying \z\ -C N, and <^ ext = <p ex t mod 2tt. 
The minima of U are surrounded by large energy barriers of 
height > 2£j , except along the special direction defined by 
a simultaneous and equal change in all variables, i.e. 9 m — 
(j)/N for all array variables. Such collective dynamics is asso- 
ciated with the black-sheep variable and has a barrier height 
of only 2_Ej . In the quantum regime, fluctuations occur pri- 
marily along this direction and motivate the use of <j> as an 
essential collective variable. 

For the Symmetric Fluxonium Model, this collective mo- 
tion of all array variables forms a normal mode for harmonic 
fluctuations around the global minimum which, for zero mag- 
netic flux, is located at <f> = 0. Anticipating the role of this 
mode we refer to it as the superinductance mode. Further 
analysis shows that the remaining N — 1 normal modes are 
degenerate and, so as to remain orthogonal to the superin- 
ductance mode, their amplitudes sum to zero mode by mode. 
We therefore call them difference modes and introduce £ M 
(/i = 1, . . . , A— 1) for their variables (Pig. |2j. The transforma- 
tion to the new set of variables {<fi, £i, . . . , £at-i} is facilitated 
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by 
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and inversely, <f> = J2 m °rn and ^ = Y^ r 
less otherwise specified sums over Greek-alphabet indices run 
over the difference modes n = 1, . . . , N — 1. Here, W is a 
[N -1)xjV semi-orthogonal matrix with components sum- 
ming up to zero in each row, J2 m W^mWum = 5^ v and 
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differs from the choice in RefP^ and proves particularly con- 
venient for the subsequent discussion of corrections from 
ground capacitances [see Eqs. |7) and ( [T0] >] which break S^v 
symmetry. After this variable transformation, circuit quanti- 
zation yields the Symmetric Fluxonium Model Hamiltonian 

#sfm = - ^E h c 'dl - 4EJ Y, d l ~ E ' cos ^ + 



(5) 



where ££' 



E^/[1 + E^/(E^N)} is equal to the black-sheep 
up to a small 1/iV-correction. 

The structure of i/sFM illustrates the utility of the 
collective-mode description: coupling between different 
modes is limited to potential energy terms. Further, at the 
relevant potential minima, all difference mode variables van- 
ish, £ p = 0, and the arguments of the array cosines [last line 
of equation |5]l] are of order 1/N. Hence, a Taylor expansion 
for small arguments can be expected to capture the essential 
low-energy physics. Keeping terms up to second order in this 
expansion obtains 
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where = (^ M /A^ — A^c^) / V2 is the ladder operator 
creating an excitation in the /i-th difference mode, A^ = 
is the oscillator length, and fl = ^/8E[Bj the 
array junction plasma frequency. 

The first three terms in the expression for Hq reproduce 
the superinductance model that was successfully used in 
Refs. [8] and |38] It describes the superinductance mode as 
the coupled system of the black-sheep junction with capaci- 
tive energy E* = E^' and a large superinductance^ L s = 
N($o/2ir) 2 /Ej with correspondingly small inductive energy 
El = Ej/N. As the second crucial insight from Hq we note 
that, within the harmonic approximation, the symmetry of the 
circuit has been extended to include arbitrary unitary trans- 
formations of the N— 1 degenerate difference modes. As the 
superinductance mode is a scalar under the action of the group 
U(A L -1) it must completely decouple from the difference modes 
in the harmonic limit. This decoupling explains, in part, the 



success of the superinductance model in matching experimen- 
tal spectra in spite of the presence of the large number of addi- 
tional degrees of freedom. The concept of symmetry-induced 
decoupling carries over to more complicated circuits that in- 
clude linear arrays of Josephson junctions. 



III. WEAK S N SYMMETRY BREAKING 

The models discussed so far obey Sjy symmetry and, at the 
level of Hq, also U(iV— 1) symmetry in the difference-mode 
subspace. To go beyond the superinductance model and pre- 
dict corrections arising from the weak interaction between the 
superinductance mode and the difference modes, we next con- 
sider mechanisms leading to symmetry breaking. We consider 
the following three mechanisms which are likely the dominant 
ones in present experimental samples: anharmonicities of the 
potential energy neglected in the above expansion (8Hjj), dis- 
order in the Josephson energies of individual array junctions 
(8Hj), and additional stray capacitances of each superconduct- 
ing island to ground (5Hc). We first derive the Hamiltonian 
expressions for each of these corrections, and discuss their ef- 
fects on the energy spectrum and eigenstates subsequently. 

We start with 5Hjj, the corrections from anharmonicities 
exhibited by the periodic potential but neglected in the har- 
monic approximation employed in Hq. Considering higher- 
order terms in the Taylor expansion of Hsfm — Hq, we find 
that the leading anharmonic corrections are given by 
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It is easy to verify that 5Hjj breaks the U(N— 1) symmetry but 
preserves the permutation symmetry under Sat. 

To derive an expression for SHj, we capture disorder in the 
Josephson energies of the array by defining Ej. m — Ej+SEj m . 
Such disorder is expected to be caused by slight variations in 
junction size and thickness, and may also be affected by junc- 
tion aging. In the absence of experimental statistics for flux- 
onium junction parameters, we choose random 5Ej m from a 
Gaussian distribution of width SEj = 150MHz and, without 
loss of generality, impose J2 m ^Jm = 0. The disorder mod- 
ifies the potential energy of the Hamiltonian, and by Taylor 
expanding we obtain 

SHj = lYl SEsm i^l N + £m W ^ 



Disorder in individual array junction parameters generally 
leads to weak breaking of both U(A^ — 1) and Sat symmetry. 

To capture corrections from stray capacitances of the su- 
perconducting islands to ground, we include the terms due to 
ground capacitances shown in equation ([TJ. Ground capac- 
itances contribute kinetic energy terms which are easily ex- 
pressed as Tj = ^(&o/2' 7r ) 2 C g j(p'j when using node variables 
ifj for each superconducting island. Assuming overall charge 
neutrality of the circuit, we can recast these additional contri- 
butions in terms of the junction variables 9 m . Accounting for 
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FIG. 3: Difference mode spectrum for 
total excitation numbers T = 0, 1, 2. 
(a) Within the harmonic approxima- 
tion, Sjv and U(2V— 1) symmetry hold 
and produce degenerate subspaces of 
dimension D T . (b) Anharmonicity 
5Hu breaks U(iV-l) but leaves S N 
symmetry intact. The irreducible rep- 
resentations of Sjv, labeled by parti- 
tions (A) or Young diagrams, give rise 
to D (A) -dimensional degenerate sub- 
spaces, (c) Corrections from disor- 
der in array junction _Bjs and ground 
capacitances, SHc + SH; break Sjv 
symmetry. The degeneracy lifting is 
shown for the T = 1 subspace. 



the ground capacitances of the two large superconducting is- 
lands surrounding the black sheep and those of the remaining 
small islands by C*, C| < Cf, C) (Fig. |TJ, the perturbation 
can be expressed as 
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where 9j =8^ and <3t — (a M — aJ I )/(v / 2A^) for [i > 1. The 
detailed derivation of equation |7]i and analytical expressions 
for the entries of the matrices M and G are provided in the 
Supplementary Material. 

After removing irrelevant global energy shifts, the effects 
of the perturbations SHjj, SHc and 6Hj can be organized into 
three categories according to their action on the superinduc- 
tance and difference modes. Perturbations may introduce cou- 
pling among superinductance states (SHI), coupling among 
difference mode states (SH^), as well as coupling between the 
two subsystems (SHl d ). 

We first discuss corrections in the 5HI category. The sim- 
plest contributions of this type are terms whose structure 
matches terms already present in Hq, and hence merely renor- 
malize the superinductance model parameters. Both SHjj and 
SHc contain corrections of this type and yield renormalized 
parameters 
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We emphasize that the inclusion of ground capacitances of 
the large islands alone leave the Sjv symmetry unharmed and 
can be completely accounted for in our model by the above 
renormalization, even in the limit that Cg is large. The only 
contribution of type SHI which goes beyond renormalization 
is the term — Ef cf) 4 /(41N 3 ) generated by SHjj. For realistic 
parameters, we find that this perturbation generates level shifts 
well below 100MHz (Fig. |4ji). 

We next turn to corrections in the SH, category with con- 
tributions from each SHjj, SHc and SHj. Perturbations like 



SHjj which do not break the Sjy symmetry, do not lift degen- 
eracies completely and necessitate the use of degenerate per- 
turbation theory. Each remaining degenerate subspace is then 
associated with an irreducible representation of the symmet- 
ric group. Our construction of the relevant irreducible sub- 
spaces can be sketched as follows. We start by decompos- 
ing the difference-mode Hilbert space into orthogonal sub- 
spaces V T with fixed excitation number T = ^ a^a^, i.e., 
Hd = Vrj © Vi £B • • ■ . In general, each V T may still be re- 
ducible under SV and should be decomposed further. 

In the decomposition, the integer partitions of N serve as la- 
bels for the irreducible representations of the symmetric group 
Sn- Here, a partition (A) = (Ax, A2, ■ ■ • , Xf) is a sequence of 
non-increasing positive integers Ai>A2>--->Af>0 
that sum to N, and are represented by a unique Young dia- 
gram. Since the inductive decomposition^ is not very practi- 
cal for 1, we instead achieve decomposition of V T by us- 
ing a restricted set of semi-standard Young tableaux 40 . (Tech- 
nical details of this procedure are provided in the Supplemen- 
tary Material.) For the low-energy part of the spectrum probed 
by experiments, we find that the excitation number T and par- 
tition (A) are sufficient to specify the relevant irreducible sub- 
spaces. 

The T = subspace, spanned by only the ground state 
|0), immediately forms an irreducible representation. The 
state is effectively bosonic and is indexed by the partition 
(A) = (N). For T = 1, likewise, Vi already forms an irre- 
ducible N— 1-dimensional subspace corresponding to the rep- 
resentation with partition (A) = (N—l, 1) (Fig. [3]). ForT = 2 
the decomposition is more interesting and results in three irre- 
ducible subspaces indexed by (N), (N—l, 1) and (N-2,2). 
The subspace labeled by (N), for example, is comprised of a 
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TABLE I: Energy corrections 
for difference mode states due 
to anharmonicities (SHu)- T 
is the total number of excita- 
tions in the difference modes 
and (A) the partition labeling 
the irreducible subspace. 
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(a) H + &H V {S N symmetric) (b) H + SHu + SH C + 5Hj (Sjv broken) 




FIG. 4: Spectrum from numerical diagonalization including effects of (a) anharmonicity, and (b) anharmonicity, ground capacitances and 
junction disorder. Dashed curves show the pure superinductance spectrum for renormalized E^' and El [equation ([8}]. Thick curves in (a) 
indicate N— 1-fold degenerate levels that remain decoupled from T = states under SHu- (b) Corrections from ground capacitances SHc 
and junction disorder SHi break Sjv symmetry. In the T = 1 manifold, ground capacitances split off several levels; smaller shifts are due to 
junction disorder 5Ef„, here chosen from a Gaussian distribution. The panels in (c), (d) show magnifications of regions marked in (b). The 
avoided crossing in (c) is primarily generated by ground capacitances. The even smaller splitting in (d) is purely generated by array junction 
disorder. Parameters are chosen consistent with experimental device (Supplementary Material). 



single Sn invariant state given by 

|T(A))^|2(^))= ^_i_ ^(at) 2 |0), (9) 

which is independent of the specific choice of W in equation 
Q. By employing perturbation theory for SH d f in each irre- 
ducible subspace V T m, we obtain the first order energy cor- 
rections SEfj[T, (A)] in Table H The resulting level shifts are 
shown schematically in Fig. [3 1. 

We next discuss the effects of ground capacitances, SH^, 
and junction disorder, SHf, on states in the in the 1 -excitation 
subspace Vi. If the ground capacitance of the big islands 
is sufficiently large compared to that of the small islands, 
E*/E* < l/N 2 , then SH d in equation is approximately 
diagonal. This decoupling is the primary motivation for our 
choice of difference modes ( Eq. [3] ). To leading order, the 
resulting energy shifts are given by 



-6E*/fi 2 {n "C N), 



(10) 



where SE d 



N 2 El 2 /(ir 2 A<:El). Similar effects from 
ground capacitances have been predicted and observed in 
Ref. [16] For higher values of //, shifts induced by Joseph- 
son energy disorder SHf become dominant in our model. For 
Gaussian distributed junction parameters SEj m , the first-order 
energy shifts SEf also follow a Gaussian distribution with 
width SEf = 6E s A 2 /2 (Figs. [5]; and^). 

Corrections in the third and final category, SH^, arise from 
coupling between superinductance and difference modes. An- 
harmonicity captured by SHff preserves Sn symmetry and 
hence, by Schur's lemma, cannot couple states belonging to 
different irreducible representations. More specifically, states 
of the form \t') ® aJjO) cannot couple to states of the form 



\£) <g> |0) under Sat symmetry, even when such states are de- 
generate (Fig. |4^). This symmetry-enforced lack of coupling 
between the superinductance mode and the lowest difference- 
mode excitations constitutes the second central result of our 
work. It is a crucial ingredient in preserving the respective 
identity of these collective modes and explains the quantitative 
accuracy of the superinductance model at low energies. The 
only difference-mode excitations that may couple to states of 
the form |^)<E>|0) are those that are bosonic, i.e., are indexed by 
the partition (N). The candidate states with lowest energies 
are <E) \(N)2) but are already well beyond the frequency 
range probed by spectroscopy in previous fluxonium experi- 
ments. 

Ground capacitances, as described by SHq, break Sjv sym- 
metry but preserve "circuit parity": mirroring of the circuit 
variables according to (pj — > tpN-j, Vext — > —fext leaves the 
Hamiltonian invariant. The superinductance mode and differ- 
ence modes with odd index (/i + 1) are odd under parity; dif- 
ference modes with even index (fi+ 1) are even under parity. 
As a result, SH^ can only couple the superinductance mode 
to every other difference mode. We find that the coupling is 
largest for small values of fi and takes the form 
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where SEg = 8E^E^N 3 / 2 /(TT 2 E h s ). Finally, all symmetries 
are broken for array junction disorder, and the resulting per- 
turbation is given by 



SHfn^SEf^ + al) 



(12) 



with SEj d ^ following a Gaussian distribution of width 
5E } At/(V2N). As shown in Figs.^d, the coupling between 
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superinductance and the // = 1 difference mode, induced by 
SHj d only, is considerably smaller than the coupling to the 
H = 2 difference mode which is dominated by SH^. 

IV. DISCUSSION AND SUMMARY 

Not surprisingly, the low-energy spectrum of the full fluxo- 
nium circuit includes, in addition to the energy levels pre- 
dicted by the superinductance model a large number of nearly 
degenerate excitations. We have identified the nature of these 
collective excitations with the difference modes at energies 
near the array junction plasma frequency. Degeneracies are 
expected to be lifted, first by U(A r — 1) symmetry breaking due 
to anharmonicity and further by symmetry breaking due 
to array junction disorder ground capacitances. The impor- 
tant consequences of these corrections include separation of 
previously degenerate levels into closely spaced multiplets. 

Josephson junction arrays provide an interesting example 
of a quantum system with many identical but distinguish- 
able degrees of freedom, resulting in representations of the 
symmetric group not readily observed in nature with indis- 
tinguishable particles. Invariance under permutations of the 
junction variables is a generic symmetry expected to be im- 
portant for any large superconducting circuit containing one 
or several Josephson junction arrays. The decomposition of 
the symmetric group Sn into irreducible representations rele- 
vant at low energies thus becomes an important tool in circuit 
analysis. For the example of the fluxonium device, we have 
shown that such symmetry strongly restricts the possible cou- 



pling between the superinductance mode, as observed in the 
experiment 8 9 , and the additional difference modes. Our re- 
sults explain the remarkable accuracy of the effective superin- 
ductance model as long as the renormalizations of El and Eq 
are taken into account, and are consistent with previous fits of 
experimental data where El and Eq were used as a fit param- 
eters, producing excellent agreemen t 8 ! 9 ! 41 ! 

The power of symmetry-based approaches in the analysis of 
future circuits is easily illustrated for the example of the flux- 
onium device. Specifically, the number of difference-mode 
states with excitations up to some threshold T grows rapidly 
as j^Ziy ■ The number of states with proper bosonic sym- 
metry, however, is dramatically smaller: For N = 43 and 
T = 5 there are 10 6 difference mode states states but merely 
6 of them possess bosonic symmetry. Harnessing exact and 
approximate symmetries of Hamiltonians for larger circuits 
will likely be a crucial ingredient in future research explor- 
ing quantum coherence in superconducting circuit networks 
of increasing complexity. 
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Supplementary Material 

Appendix A: Irreducible representations for difference modes 

In this appendix we discuss the decomposition of the 
difference-mode Hilbert spaces V T into subspaces that trans- 
form irreducibly under Sat symmetry. Since some of the math- 
ematical tools employed may not belong to the physicist's or- 
dinary repertoire, we provide definitions along with concrete 
examples where appropriate. In terminology and notation, our 
discussion closely follows, the excellent book by SagarRH 

The subspace V T comprises all difference-mode states with 
total excitation number T. It is spanned by the orthogonal 
states 

where we assume weakly ordered mode indices 

/ii<---<MTG{l,2,--- ,N-1} 

to avoid double counting. The Sn symmetry displayed in the 
ideal fluxonium circuit pertains to permutations a € Sat of 
the array junction variables 61,62, ■■■ ,6n- Such permuta- 
tions also induce linear transformations in the operator space 
C{a{, a\, . . . , a] v _ 1 } spanned by the difference-mode cre- 
ation operators. To understand how at transforms under per- 
mutations, we recall the definition of the creation operators in 
terms of junction variables: 

a l = WftmiOm/At - A^dg m )/V2. (A2) 

m 

Using the identity ^ W^W^n = S mn - i, one finds that 
the difference-mode creation operators transform according to 

<r(al) = Y,S(*) lu ,al. (A3) 

V 

We remind the reader that, by our convention, sums with Latin 
(Greek) summation indices always range from 1 to A^ (1 to 
A" — 1). The (AT — 1) x (AT — 1) transformation matrices are 
given by 

J2 W ^mWu„D(a) mn . (A4) 

Here, D(a) mn = S miCr / n ) denotes the N x N permutation 
matrix for the group element a £ Sjy. (The matrices D(a) 
form the defining representation of Sat.) 

By the relation Vi = C{a{, a 2 , . . . , ajy_x}|0), the transfor- 
mation matrices S(cr) in equation ( |A3[ ) define an orthogonal 
Sjy representation of degree A" — 1 in the one-excitation sub- 
space. Similarly, for higher excitation numbers T > 1, the 
group action for products of creation operators, 

(A5) 



determines the representation of Sat in the subspace V T . Given 
these representations, our central task is to decompose each 
V T into its irreducible subspaces. As an aside we note that 
in the special case of Vi, simple arguments based on group 
characters can be used to show that Vi is already irreducible 
and coincides with the irreducible representation indexed by 
the partition (AT— 1, 1), for which group characters are known 
to be txD(a) - 1 for arbitrary N (Ref.gO| section 2.12). 

The common approach for decomposition of such product 
representations is inductive and requires successive decompo- 
sitions for Si, S2, ■ ■ ■ , Sat, see, e.g., Ref. [39] For large N, 
however, that strategy is not very practical. Following the 
treatment by SagarPa we thus employ an alternate approach 
using a restricted class of semi-standard tableaux. (We explain 
the meaning of these words in due course). 

As our first step in constructing the decomposition of each 
V T , we define the pseudo-creators b\ for n = 1, 2, . . . , N by 

b{ = J2 W ^ a l (A6) 

= [(6 n - &/N) /A s + A £ {d 6n - d+)] /V2. 

As one would expect, pseudo-creators w n increase the total ex- 
citation number in the difference-mode subspace by one. The 
number of pseudo-creators, however, is A^ and thus exceeds 
the number of difference-mode creation operators by one. 
Indeed, the pseudo-creators obey J2 n = an d, hence, are 
not linearly dependent. They obey the non-standard commu- 
tation relation 

[b m ,bi] = 5 mn -l/N. (A7) 

For the price of this anomalous commutator, we obtain op- 
erators which transform with elegant simplicity. Specifically, 
under array variable permutations a 6 Sat, the w n operators 
simply undergo permutations: 

a{bi)=Y J D{a) nm b\ n = bl {n) . (A8) 

This simple transformation law will be crucial for finding the 
irreducible subspaces of V T . 

We next extend the language of difference-mode excitations 
to pseudo-mode excitations and define the states 

N 

|t) = Y[ (O*" 1°) (not normalized), (A9) 

n=i 

where the vector t = (t\, ti, ■ ■ ■ , t^) specifies the excitation 
numbers t n € No for each pseudo-mode b^. Using the inverse 
of equation ( |A6| l, = J2 n W^nb'n, it is simple to confirm 
that 

V T = span{| t ) I ^ n t n = t}, (A10) 

i.e., the pseudo-mode excitations span the difference mode 
subspaces V T one by one. 

For a state with given pseudo-mode excitations t = 
(ti,t2,- ■ ■ ,tN), equation ( |A8| l implies that a E Sn simply 
permutes the pseudo-mode excitation numbers according to 

t -t t' = (* CT (l),i<7(2)j • • ' ,t<r(N))- 
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(a) mode: n=l 2 3 N-l N 



ti = 2 t 2 = 1 f 3 = t N .i= 3 t N = 1 

(b) exc. #: T T-l 2 1 




FIG. 5: Labeling schemes for difference mode excitations, (a) La- 
bel specifies excitation numbers, mode by mode, (b) Label specifies 
mode indices, excitation number by excitation number. 



Example: The state with 3 excitations in pseudo- 
mode n = 1 and 4 excitations in pseudo-mode 
n = 2 is t = (3, 4, 0, . . . , 0). A permutation 
may transform it into t' = (4, 0, 3, 0, ... , 0), for 
example but not into (2,5,0,.. .,0), even though 
the latter state still has the same total excitation 
number T = 7. 

For each state |t) with T = J2 n we define the subspace 
spanned by itself and its permuted partner states: 

V [t] := span{|t') = \at) \ a e S N } C V T , (All) 

As suggested by the notation, [t] may be understood as an 
equivalence class when defining t ~ t' :<4> there exists a a £ 
Sn such that t' = at. By construction, the Vrt] form Sn- 
invariant subspaces, and their unions cover each V T 

V T = |J V [t] . (A12) 

{M|T=£„t„} 

Note that, due to linear dependence of the pseudo-modes, sub- 
spaces for inequivalent excitation classes [t] n [f] = may, 
nonetheless, have a non-zero intersection, V[t] H V[f] 7^ 0- 

We will first discuss the decomposition of V[t] into irre- 
ducible subspaces as if b n were orthogonal modes. In step 
1, we thus drop the cautionary prefix "pseudo" temporarily 
and show that the basis vectors 1 1 ) can then be relabeled in 
such a way to reveal isomorphism between Vrt] an d me cor ~ 
responding permutation module M ' . In step 2, we then uti- 
lize the important theorem for the decomposition of M Xt that 
identifies semi-standard tableaux as the indexing set for all 
irreducible subspaces. In both steps, we introduce the nec- 
essary terminology and explain the construction. We do not 
provide proofs of the underlying theorems but refer the inter- 
ested reader to Sagan's boolf^, chapter 2. Finally, in step 3 
we return to the issue of linear dependence of b n modes and 
show how the usual construction can be modified to account 
for the linear dependence in a simple fashion. 



STEP 1: Isomorphism between V[ t ] and the permutation 
module M At 

The excitation numbers t = (ti, t%, . . . , tjy) label each 
state in V[ t ] by specifying excitation numbers, mode by mode. 
An alternative labeling scheme (Fig. [5]) consists of specifying 
mode numbers, excitation level by excitation level. 

To specify the procedure of switching from the first labeling 
scheme to the second, we define the partition (At) associated 
with t as follows. Consulting t and for each integer / = 
0, 1, . . . , T, count 

X'f = (number of modes with / excitations). (A13) 

The resulting sequence (X' , X[, . . . , X' T ) sums to N, the to- 
tal number of modes. By sorting entries in this sequence in 
decreasing order and dropping all entries, we obtain the par- 
tition (A t ) associated with t. Excitation numbers t,t' in the 
same equivalence class always have the same associated par- 
tition. 

Example: for N — 4 and excitation num- 
bers t = (1,3,0,1), one obtains the sequence 
(1, 2, 0, 1, 0, 0) and thus the associated partition 
(At) -(2, 1,1). 

The partition (At) = (Ai, A2, . . . , Xf) is represented by 
a Young diagram: an array of squares where row / has A/ 
squares. From the Young diagram (At) we obtain the Young 
tableau 9 t associated with t by filling the boxes with the 
mode indices from 1 to N in such a way that mode indices 
with the same excitation number appear in the same row. 

Example: Fort = (1, 3, 0, 1), the Young diagram 
of (At) = (2, 1, 1) and a Young tableau associ- 
ated with t are given by 



(At) = 







1 


4 




and O t = 


3 








2 





From the sorting function used to order the entries in the 
partition (At), one can infer which row in the tableau refers to 
which excitation number. As a result, the Young tableau lists 
the mode indices corresponding to each excitation number as 
intended. One additional modification is required to turn it 
into the desired state label. 

For this, note that the transposition (interchange) of two 
mode indices with the same excitation number leads to a new 
tableau (consider interchange of the entries 1 and 4 in our ex- 
ample above) but not to a new state 1 1 ) . To remove this ambi- 

guity, we define two tableaux as row-equivalent, 0t ~ 8^ :<^» 
permutations of elements within each row can transform Q' t to 
0t. The resulting equivalence class [0t] is called the tabloid 
associated with t and serves as the new label for each state. 

Example: Using t = (1, 1, 0, 3) as above, the as- 
sociated tabloid is 
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It is useful to note that [0 t ] can also be expressed as [0t] = 
{<70t | u £ R0t} where R0t C Sn is the subset of permuta- 
tions which only interchange entries in each row of the tableau 
©t. R0t called the row-stabilizer of t . Below, we will also 
encounter the column-stabilizer C0 t with the analogous defi- 
nition referring to columns instead of rows. 

With this, we have established a one-to-one map 1 1 ) O 
I [0t] ) which achieves the important goal of relating V[ t ] to a 
central object in the representation theory of Sn- the permu- 
tation module M At defined by 



M At = C { [0t] 1 1 G [t] } . 



(A14) 



Since each permutation module is defined in terms of tabloids, 
the group action for V[ t ] and for M * is easily verified to be 
identical, and the two vector spaces are hence isomorphic as 
Sjv representations. 



STEP 2: Decomposing V[t] and constructing basis vectors for all 
irreducible subspaces 



With this preparation, one now obtains the basis states span- 
ning the instance(s) of within Vrt] from 



1ST; ©ft) 

= E E sga(a)ar b\ 



f 0(1) 1 0(2) 
& 2 



t #W 



N 



(A16) 

10) 



where permutations a, r £ Sjv act on the mode indices, i.e., 
the subscripts of the b^ n operators as before. (Caveat: as de- 
fined above, the basis states are not normalized yet.) Each 
semi-standard tableau £^ yields an irreducible subspace 



span 0£) | all standard tableaux 0^} . (A17) 



By considering all possible partitions (/i) and associated semi- 
standard tableaux , we thus completely decompose V[t] 
into linearly independent irreducible subspaces. (The set of 
partitions one needs to consider can be restricted by consider- 
ing dominance ordering of partitions^.) 



The great benefit of identifying V[t] as isomorphic to M Xt 
lies in the availability of mathematical tools for decomposing 
the permutation modules into their irreducible subspaces (see 
Sagan 40 , section 2.10). For the symmetric group Sat, each ir- 
reducible representation 5 M is labeled uniquely by a partition 
(fi) of N. Consequently, the general decomposition takes the 
form 



M- 



(A15) 



where £ No is the multiplicity of the irreducible sub- 
space within M x . For a given Vm — M At , we wish to 
obtain the basis vectors for each of the copies (if any) of 5 M 
contained in it. The basis vectors are obtained by means of 
semi-standard tableaux, which we define next. 

The associated semi-standard tableau is constructed 
from the Young diagram for (/i) by filling its squares with the 
excitation numbers t±, t2, ■ ■ ■ , tjv in such a manner that entries 
in each row weakly increase (<), and entries in each column 
strictly increase (<). Along with ££", we consider standard 
tableaux 0g of the same shape, which is a Young tableau with 
entries increasing in each row and column. This way, we can 
set up a lookup function d : {1,2,..., N} — > {ti, ta> • • • jijv} 
that extracts the position of the integer n in the reference 
tableau 0g and returns the excitation number entry found in 
the semi-standard tableau at the corresponding position. 
To keep notation simple, we usually do not make the depen- 
dence of d on and 0g explicit. 

Example: For [t] = [(0,1,0,3)], the associated 
semi-standard tableau with shape (/z) = (2, 2) 
and the accompanying reference tableau are: 









1 


3 



M 



1 


3 


2 


4 



The lookup function then yields the results: 

12 3 4 



STEP 3: Linear dependence of b n modes and restricted 
semi-standard tableaux 



In steps 1 and 2 we have ignored the linear dependence of 
b n expressed by the constraint b n = 0. Once linear de- 
pendence is taken into account, the states from Eq. ( |A16| > still 
span, 



V T = span 



££; eg) 



0g semi-standard and 



standard tableaux, T = J2 n ^ 

(A18) 

However, they are linearly dependent in general. Fortunately, 
removing this linear dependence can be achieved by a harm- 
less modification of our previous procedure. This modifica- 
tion consists of an additional restriction on the set of admis- 
sible semi-standard tableaux Specifically, we will show 
that restricting the semi-standard tableaux to the set 



K T = | has no "1" in its first row, J2 n l n = T }, 

(A19) 

removes the linear dependence and 

Bj = 9£) | £ TZ T , 0£ standard tableau} (A20) 

forms a basis for V T . Furthermore, each subspace spanned by 
states with a fixed restricted semi-standard tableau retains 
its character as an irreducible representation indexed by the 
partition (//). 

To prove this assertion, we first show that every state 
l^t'i eft) obtained for a "forbidden" semi-standard-tableau 
££ ' ^ TZt can be written as a linear combination of states 
I > ef) from the restricted set, i.e. E^ 1 £ 1Z T . Consider 
the vectors constructed in Eq. ( |A16| l and note that the row- 
stabilizer can always be separated into the stabilizer of only 
the first row Ri and the stabilizer of all remaining rows R': 



10 3' 



TGRei' t'GR' tGRi 



(A21) 
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Example: In this and all following examples we 
consider the "forbidden" state vector 






1 


2|. 


1 


3 


5| 


1 


4 




2 


4 





The stabilizer for row 1 consists of Ri = 
{e, (13), (15), (35), (135), (153)}; the stabilizer 
for the remaining rows is R' = {e, (24)}. 

Proceeding with the decomposition of the "forbidden" state 
vector |E{* ; O^) in terms of states with restricted semi- 
standard tableaux, let M = {mi,m2, ■ ■ ■ ,tojv} denote the 
entries of the standard tableau (reading left to right, row 
by row), and M\ = {mi,m2, • • • , to,.} the entries in row 1 
only. Similarly, let (ii^r • " >t r ) be the integer excitation 
numbers in the first row of the semi-standard tableau E£ and 
(without loss of generality) assume that t\ = • • • = t q -i = 
and t q = 1 for q < r. 

Next, we introduce the sets A m = {n r , n r -i, ■ ■ ■ , n m } and 
rewrite the stabilizer of row 1 as 

E T K^ti-eJ (A23) 

= (?-i)« Ee E fc--- E ^- 

n r GA/i n r _i£Afi\A r n q £ Ah \A g + 1 

By construction, each successive sum over pseudo-mode in- 
dices n r , n r _i, ... is associated with weakly decreasing ex- 
citation numbers t r > £ r _j > ■ ■ ■ and the final sum over 
pseudo-mode indices n q corresponds to case of a single exci- 
tation (entry "1" in row 1). 

Example: Continuing with our previous example 
the above equality takes the form 

E-( 6t ^ 2 )= E € E b l 

i-eRi n 3 e{l,3,5} n 2 6{l,3,5}\{n 3 } 

Next, we use the linear dependence of the pseudo-modes to 
rewrite the final sum in Eq. ( |A23) > as 

E b k = - E 6 "- E b l (A24) 

n q £M 1 \Ag +1 n£A q + 1 neM\Ah 

The transformed expression has two separate sums over n: 
a sum over pseudo-modes that, according to Eq. ( |A23| >, are 
already occupied, and a sum over pseudo-modes in rows 
2,3,... of the standard tableau. The increase of excitation 
number produced by Eq. ( |A24[ ) hence only affects pseudo- 
modes that are already occupied in Eq. ( |A23[ ). 

Example: Again continuing with our previous ex- 
ample we find 

E b l = - J2 b n- E & t < A25 > 

nae{l,3,5}\{n 3 } ne{n 3 } nS{2,4} 



Finally, we inspect the full state vector by using Eqs. ( | A2 1] >, 
( |A23[ ), and ( |A24[ ). The terms resulting associated with a sin- 
gle index n [Eq. ( |A24| >] can be re-expressed as a sum over 
the complete row stabilizer R and associated with a tableaux 
where the entry "1" in the first row of E£* has been eliminated 
and another entry 1 < t n of E^ 1 ' has been increased by one: in 
such a way the state with "forbidden" semi-standard tableau 











1 


tq+1 






tr+1 




t r2 











is decomposed into a linear combination of states with 
tableaux 














tq+1+1 




tr 


tr+1 




t r2 





















tq+1 




t r + l 


tt+1 




t r2 





















tq+1 




tr 


tr+1+1 




t r2 





















tq+1 




tr 


tt+1 




tr 2 +l 











In cases where the above procedure results in a tableau that is 
not semi-standard, a straightening algorithm can be applied to 
generate the corresponding semi-standard tableau 40 , section 
2.6. Importantly the straightening algorithm does not change 
the content of the tableaux and thus, does not change the fact 
that our procedure expresses the "forbidden" tableau in terms 
of semi-standard tableau with fewer "1" entries in row 1. Us- 
ing this procedure, repeatedly if necessary, we can decompose 
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any state with "forbidden" semi-standard tableau as claimed. 
It is important, of course, that the removal of states occurs at 
the level of entire subspaces (indexed by the forbidden semi- 
standard tableaux) while the group action and hence the irre- 
ducibility of the remaining subspaces is unharmed. 

Example: We complete our running example by 
decomposing the "forbidden" state vector 






1 


2|. 


1 


3 


5| 


1 


4 




2 


4 





(A26) 



in the restricted basis. Following our previous 
steps, the state vector can be expressed as 



E E s S n ( cr ) (JT 
crec tgR' 



(A27) 



- E C - E < E * 

n 3 eMi n 3 £Mi rie{2,4} 

Together with a combinatorial factor (here: 1 /2), 
the sums over 713 and r' G R' can be recombined 
into the full row stabilizer: 



-jEE sgn ( CT ) <jr 



ne{2,4} 



We have thus completed our goal of expressing 
the original "forbidden" state vector as a linear 
combination of the "restricted" state vectors 









3|. 


1 


3 


5 


1 


4 


1 


2 


4 











2. 


1 


3 


5 


2 


4 


1 


2 


1 











2 . 


1 


3 


5 


1 


5 




2 


1 





With the decomposition of "forbidden" semi-standard 
tableaux in hand, we conclude by showing that the states 
l^t ! ®st) w i m restricted £^ 6 71% not only span each V T but 
are linearly independent. The proof is based on a simple di- 
mensional argument. Counting the number of basis elements 
with T excitations, we find 



dim V T 



(JV-2 + T)! 
t!(AT-2)! 



(A28) 



This should be compared with the dimensionality of sub- 
spaces constructed with restricted semi-standard tableaux 
€ TZ T . For ^ n t n = T < 5, the explicit listing of re- 
stricted semi-standard tableaux is given in Table II. For com- 
parison with Eq. ( |A28| >, note that the dimension of the ir- 
reducible representation of Sn indexed by partition (/x) can 
be obtained by the hook length formula, see Sagan 40 , section 
3. 10. We have verified that 



dp = 



(AT-2 + t)! 
t\(N- 2)! 



(A29) 



for all T < 5, and leave it as a conjecture that equality contin- 
ues to hold for all higher T. 



Example: As an example we consider T = 2. 
Using the hook-length formula the sum of the di- 
mension of the irreducible subspaces indexed by 
the semi-standard tableaux 



I I - I I 2 I , 






- 1 | 










2 




1 


1 





(A30) 



is equal to 1 + N - 1 + N{N - 3)/2. Sim- 
ple arithmetic shows this is equal to dim V T= 2 = 
N(N - l)/2 as expected. 

In summary, the subspaces indexed by restricted semi- 
standard tableaux decompose each V T into its irreducible sub- 
spaces. As a final remark, we note that for T > 2 multiplicities 
of irreducible representations can exceed 1. (Through T = 5, 
the largest multiplicity that occurs is 5, see Table II.) In such 
cases, the usual Gram-Schmidt procedure may be employed 
to generate orthogonal irreducible subspaces. 



Appendix B: Calculation of perturbative shifts according to 
irreducible subspaces 

We next discuss the calculation of first-order shifts of 
energy levels in each irreducible subspace under the Sn- 
symmetric perturbation SHjj. In general, irreducible sub- 
spaces for the difference modes are labeled by restricted semi- 
standard tableaux. (For T < 2, is sufficient to specify T and 
the partition (A) instead of full-blown semi-standard tableaux. 
In this appendix, we continue to employ the restricted semi- 
standard tableau notation.) 

To calculate the first order shifts, we choose a unique el- 
ement from each subspace. This is done by fixing a refer- 
ence standard tableau 0* f , which we choose as the standard 
Young tableau of shape (A) with entries 1 through N filled in 
column by column. Using this reference tableau, we obtain 
one representative state in each irreducible difference-mode 
subspace, which we denote by ) (£^ £ 7l£)- 

Example: The reference standard tableau for the 
partition (A) = (N 



^ref — 



2 


2) 


is 


1 


3 


5 | - \N 


2 


4 





(Bl) 



The state acting as the representative for the T = 
2, (A) = (A— 2, 2) subspace is then given by 









1 - 1 


1 


1 





£5>gn(a)<TT (44)10) 

<t6C t£R 
bl - b\) (b\ b\ 



|0>, 



where the column and row stabilizers C and R 
are defined with respect to the reference standard 
tableau. With Eq. ( |A7) , we obtain 









1 - 1 1 


1 


1 





5Hfj 









1 - 1 1 


1 


1 





2££(1 - 1/A) 



for the first-order shift of the irreducible subspace 
due to the effect of anharmonicity. 
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Table II: Irreducible subspaces for difference modes with total excitation number T < 5. In this table the equiva- 
lence class [t] with A^ pseudo-modes with exactly / excitations is denoted A °l Al ••• where all entries with A^ = 
are suppressed. All semi-standard tableau of shape (fi) and content [t] are listed in the corresponding row and column. 



T 


[t] 

A ol A i • • • 


(u) = (N) 


(iV-1,1) 


{N-2,2) 


(iV-2 l 2 ) 


(N-3 3) 


(N—3 2 l 2 ) 


(N-A 4) 


(iV-4 3 1) 
i jv ^) u ? x ^ 




n 


nJV 
u 


1 o 1 - - 1 o 1 


















1 


o N - 1 i 






o --In I 
i_ 
















2 


N - 1 2 


1 •■• 1 1 2 




n — I o I 

2_ 


















o N - 2 i 2 







1 


... o| 

1 














3 




1 1 •■• 1 1 3 




o -lol 

3_ 


















N - 2 12 










L 
I 












[0j-|0|2| 





...o| 




1 


2 















1 


> - o| 

L 1 










4 


Cr M 


ol-lok 




o --|o| 


















0^-213 










L 

i_ 












|0j-|0|3| 





... o| 




1 


3 




JV-2 2 2 







-I0I2I 


(} 


... o| 














1 1 1 1 2 1 2 1 






2 


2 




QiV-3^2 















1 
_2 


1 


ol 











-■ 0I2I 





)... o| 


1 


1 


1 


L 2 




qAT-4^4 

















1 




1 




1 


1 


ol 






5 




1 1 ■-- 1 1 5 




n --|o| 


















0^-214 










L 
i 












[oj-loUl 





...o| 




1 


4 




(^-223 







2 


... o| 

3 




a —|o| 












1 1 1 1 2 1 3 1 



3 


-I0I21 




2_ 
I 




qW-3120 

VJ _L O 















1 
_3 


1 


ol 











■■ o|s| 





j...o| 


1 


1 


1 


L 3 




W-3 12 2 




|0j-l0|2|2| 




1 


-■ 0I2I 

2 



[2 


-I0I2I 



1 


j...o| 

I 2 





1 

_2_ 


2 


ol 










W - 4 1 3 2 












1 1 


■■■ 0I2I 
1 








1 




1 




1 


■2 


ol 


- 0] 

111 

L2J 
























0-3 

11111 



14 



Appendix C: Lagrangian for superinductance and difference 
mode variables 



After transforming to superinductance and difference mode 
variables, the Lagrangian of Eq. ([TJ can be cast into the form 



SFM 



16.E 



(CI) 



This expression for the Lagrangian has important advantages 
over its equivalent form expressed in terms of m . First, the 
kinetic energy is now diagonal. Second, low-energy minima 
of the potential energy U are in locations where each differ- 
ence mode variable vanishes, £ p = 0. Third, fluctuations 
between minima are dominantly described by the <fi variable. 
The ability to simultaneously expand around £ M = for each 
local minimum of U is key in the derivation of the superin- 
ductance model used previously^E. 



Appendix D: Incorporating capacitances to ground 

The (N + 1) node variables ipj can be expressed in terms 
of the Adjunction variables 6 m when using the constraint that 
the total charge . rij of all superconducting islands be zero. 
To see this, we may use t = ipo as a reference variable and 
express every other node variable ipj (J = 1, . . . , N) as 



(Dl) 



Note that f$o/27r represents a uniform voltage shift of all su- 
perconducting islands relative to ground, and that r is cyclic, 
i.e., the Lagrangian is independent of r. Hence, its conjugate 
momentum is conserved: dC/di = const. This constant of 
motion, in fact, corresponds to the total charge n tot since 



dC 
df 



E 



difj dC 
dr difj 



^tot- 



(D2) 



Imposing the constraint n tot = thus allows us to eliminate f 
from the Lagrangian, and to work with a Lagrangian (strictly 
speaking, a Routhian) which only depends on 9, n and m . Us- 
ing this procedure the contribution to the kinetic energy due to 
ground capacitances i($o/27r) 2 Yli=o CgiVtf m Eq- (D takes 
the form \ Y.mn Qmn,0 m n where 



Gmn 



(<£ /27r) 



2 min{m,n} — 1 TV 



2^i=0 i=0 j=max{m,n} 



(D3) 



Then, assuming the two large superconducting islands sur- 
rounding the black sheep junction have ground capacitance 



Cg = e 2 /2_Eg while the ground capacitance of the remain- 
ing small array islands are Cg = e 2 /2E* and using the vari- 
ables £p defined in Eq. ([3]) the kinetic energy terms of the 
Lagrangian [Eq. ([TJ] takes the form 



16 



N-l 



{M^ + G^tniu. 



(D4) 



Here for compactness we use the shorthand £ = <f>. The 
symmetric A^x A^-matrices M and G have the following form: 

M 00 = 1/E%, Afoju = 0, 
and 

(N- l)(JV-2) 



M^ = 5^ v /E* (D5) 



Gn 



Gov — - 



1 

~2Ej 
o 



VINE?; 



(D6) 



Guv " 



2E*V2N s(/z) 2 
4S|s( M ) 2 



EgO,, o„ c(ju) c(i/) 



NEf[2 +{N- l)E b g /E*]s([i) 2 s(v) 2 

where c(/j.) and s(/j.) are shorthand for cos(itij,/2N) and 
sin(7r/i/2A^), respectively. Furthermore, the coefficient o M is 
1 whenever the index // is an odd integer, and zero otherwise. 

Performing the Legendre transform, the perturbation from 
ground capacitances takes the form 



N-l 

6H C = 4: [( M ' 



G)- 



M~ 



For small ground capacitances, the entries of G are small com- 
pared to those of M and we approximate (M + G) _1 ~ 
M^ 1 — M _1 GM _1 , which yields equation ([7]) in the main 
text. 



Appendix E: Limits of the perturbative approach 

Corrections in the main text are treated perturbatively and 
we briefly comment on necessary conditions for this approach 
to be valid. First, we remark that the energy scales 5Ej d and 
5Eq of Eqs. ( fPTj ) and ( fT2| ) must remain sufficiently small rel- 
ative to the typical energy scales of the superinductance spec- 
trum. Secondly, when the magnitude of (0\5H^\(N)2) or 
(0\5Hc(a\) 2 \0) becomes of order 2il, the ground state of the 
difference will requires non-perturbative corrections. To pre- 
vent this, the respective inequalities 



sfN]2 < 16/Af , (ATA 5 /27r) 2 < ££/J3£ (El) 

must hold. Thus to connect with the A^ — > oo limit (see, e.g., 
Ref. l38l) a different approximation scheme to model the low 
energy spectrum of fluxonium may become necessary. How- 
ever, in the case of the fluxonium samples previously studied 
in Refs. |8l9l the range of validity of the perturbative approach 
is well satisfied. 



Appendix F: Parameter values used in numerical calculations 

The specific parameters used in all calculations in the main 
text are N = 43, E' d c = 1.0, Ej = 26.4, E b c = 3.8, E* = 8.9, 
SE] = 0.15, E\ = 5 and E't = 250. Using Eqs. |g) these 

parameters yield E$l — 2.5 and El = 0.53; all energies in 
units of /iGHz. 



